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We consider fermions on a 2D square lattice with a finite-range pairing interaction, and obtain 
signatures for unconventional pair-symmetry states, (1^2 _y2 and extended-s (s*), in the BCS-BEC 
crossover region. We find that the fermion momentum distribution function, vf.^ the ratio of the 
Bogoliubov coefficients, Vk/uk^ and the Fourier transform of vf. are strikingly different for d- and s* 
symmetries in the crossover region. The chemical potential and the gaps functions for both pairing 
symmetries show several interesting features as a function of interaction. Fermionic atoms in 2D 
optical lattices may provide a way to test these signatures. We discuss current generation cold atom 
experiments that may be utilized. 



In recent years there have been fascinating discover- 
ies of several classes of systems exhibiting many-body 
states with paired fermions, that may have unconven- 
tional pairing symmetry, different from that in s-wave 
BCS case. These range from heavy fermions, high-Tc 
cuprates, iron-pnictides to ultracold fermions. There are 
intense theoretical and experimental efforts aimed at de- 
ciphering pairing symmetries or mechanisms in these sys- 
tems. 

Ultracold neutral atoms [1 present an unprecedented 
opportunity to study the physics of quantum many- 
particle systems. Subjected to positive and negative de- 
tuning using Feshbach resonance technique, these pro- 
vide realizations of weak-coupling BCS to strong coupling 
BEC crossover behavior, and the unitarity limit, where 
the scattering length is infinite. 

Cold fermionic atoms in optical lattices JT, ^ consti- 
tute an intriguing set of systems. Tunibility of filling 
(particles per site), hopping kinetic energy or interparti- 
cle interaction render optical lattices unprecedented test- 
ing grounds for models of correlated electron systems. It 
has been suggested [4] that atoms in optical lattices, con- 
fined to the lowest Bloch band, can be represented by the 
archetypal condensed matter Hamiltonian, the Hubbard 
model with hopping t between neighboring sites, and 
on-site interaction U . Interatomic magnetic or electric 
dipole interactions, and possibly multi-band couplings, 
can result in finite-range interactions, that could simu- 
late the extended Hubbard model [5j . Duan [6 has shown 
that on different sides of a broad Feshbach resonance, the 
effective Hamiltonian can be reduced to a t-J model, fa- 
miliar in studies of correlated electron systems, J being 
related to spin or magnetic coupling. Calculations 0^9] 
predict that attractive-U Hubbard model gives rise to s- 
wave superconductivity, while the repulsive-U model re- 
sults in an antiferromagnetic or a (i-wave superconduct- 
ing phase depending on filling. There has been sugges- 
tion [10 that the underlying physics of the high Tc super- 
conductors may be understood by studying optical lattice 
systems. Also, due to the possibility that the cuprates, 
possessing short coherence lengths, could fall in the BEC- 



BCS crossover region, the crossover problem [TTJ [12] has 
received considerable theoretical attention. Several au- 
thors employed continuum models [MHTT] , focussing 
mostly on conventional s-wave pair symmetry. Lattice 
models with on-site or nearest-neighbor attractions have 
also been considered [3 [HI J.7-19 . There is also a large 
body of theory work [20l [21] specific to cold fermions. 

Observation of Mott-insulator behavior [22 and s-wave 
superfluidity (for attractive interactions) [23 in 3D opti- 
cal lattices represent remarkable feats. Current searches 
for new phases, such as antiferromagnetism (AFM) or 
unconventional pairing are greatly facilitated by recent 
development of analogs to existent powerful experimen- 
tal tools in condensed matter physics; examples are mo- 
mentum resolved radio frequency (rf) and rf pairing gap 
spectroscopies |24l [25] , tomographic rf spectroscopy [27] , 
out of lattice time-of- flight measurements [26^, fluores- 
cence imaging [28 , and momentum-resolved photoemis- 
sion spectroscopy [29], analogous to angle resolved pho- 
toemission spectroscopy (ARPES) in condensed mat- 
ter [30 . 

In this paper, we study zero temperature fermion pair- 
ing in a 2D square lattice in the BEC-BCS crossover 
regime using a finite-range pairing interaction. As repre- 
sentative cases of unconventional pair symmetry, we con- 
sider two even-parity representations of the cubic group, 
namely the £ = 2 dx'2-y2-wme^ and the ^ = extended s- 
wave (s*). We present several new results, in particular, 
specific signatures of states with unconventional pairing 
gap symmetries as one goes between weak and strong 
coupling regimes. We expect our results to be of rele- 
vance to fermionic atoms in 2D optical lattices, and to 
correlated fermion model systems these would simulate. 
A key result is the remarkable behavior of the fermion 
distribution function, v^, (related to momentum distribu- 
tion, n/c): For a d-wave pairing gap function, v'^ changes 
abruptly from exhibiting a peak at the Brillouiun zone 
(BZ) center (0,0) to a vanishing central peak accompa- 
nied by a redistribution of the weight around other parts 
of the BZ ((0, ±7r),(±7r, 0)) as the system crosses from the 
weak-coupling BCS to the strong-coupling BEC regime. 
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Its Fourier transform exhibits a "checkerboard" pattern 
in real space. By contrast, v'^ changes smoothly in the 
5* -wave case. We find similar signatures in the ratio of 
Bogoliubov coefficients v^/uk^ related to the phase of the 
superffuid wavefunction. 

Our finite-range pairing interaction is obtained from 
the extended Hubbard model for two equal species pop- 
ulation system on a 2D square lattice: 

H = (-tc+c^-^ + H.c.)+U^ni^ni_^ 

<ij>a i 

- V ^ nicjUja' - fio^rii, (1) 

<ij>aa' i 

where t is the hopping, jHo the unrenormalized chemi- 
cal potential, U the on-site repulsion and V the nearest- 
neighbor attraction, a is the "spin" index, that could 
refer to hyperfine states in optical lattices. In mean-field 
theory, the Hartree self-energy terms renormalize jHo such 
that /i = /io + + l^vif) where fiu{f) and 

are filling-dependent corrections to /i. We work with the 
renormalized /i so as to properly deal with weak and 
strong couplings, and take /J^j^if) = Ji/, where Ji = /7, 
-V. The fining / = N/2M, with TV the number of par- 
ticles, M the number of lattice sites, and the spin degen- 
eracy factor of 2. In correlated electron systems, interac- 
tions are mainly Coulombic in origin, and V is typically 
an order of magnitude down from U. These are scaled 
by t for a convenient characterization of weak and strong 
coupling, and to be broadly applicable. In optical lat- 
tices, an effective Hamiltonian, similar to Eq. (1), could 
be deduced with V arising from dipolar, or multi-band 
couplings. 

On Fourier transforming and retaining interactions be- 
tween particles with equal and opposite momentum, the 
reduced pairing Hamiltonian assumes the form: 

Hpair = ^{^k - l^)clck + ^ Vkk'C'l^ctf.^C-kCk (2) 
k kk' 

where in the tight-binding approximation, = 
— 2t {cos kx-\- cos ky); Vkk' = VQ{cos{kx — k'^)^cos{ky—k'y))^ 
which is non- separable. Using the standard BCS vari- 
ational ansatz, \^bcs >= nk('^k + ^kCkC^k)!^ 
obtain the T = gap equations for the gap functions 
^d,s _ /\^(^yj(^(3os /Ca; ±cos/c^) with dx'2-y2{-) and 
symmetries, 

1 1 cos k x(cos kx ^ cos ky) 

where E^''* = ((e/e - /i)^ + A^(cos /c^ ± cos fc^^)^)^/^ The 
Bogoliubov coefficients are given by, 

Kl^ KP = ^(i±|^)- (4) 



The ratio = -(e/e-/i))/A^'^* . We readjust 

/i for strong attractions by supplementing the T = gap 
equation with the number equation [12 : 

BZ 

k ^k 

This determines the self-consistently readjusted /i, which 
is no longer ffxed at the Fermi level, and makes the 
gap equation applicable over the entire range of fflling, 
thereby the BCS and BEG regimes. To allow for strong 
scattering, sums are performed over the entire BZ. The 
natural momentum cut-off afforded by the lattice avoids 
any possible ultraviolet divergences. 

Remarkable differences stem in an essential way from 
differences in gap symmetry. The dx2-y2 gap van- 
ishes along the lines i/c^^ = in the 2D BZ, i.e. at 
four points on the Fermi surface(fs), the location of which 
depends upon filling. The 5* gap A| coincides with the 
tight-binding fs at exact 1/2-filling, and is nodeless oth- 
erwise. Here, /i < 0, with /i = —At at the bottom of the 
band. Owing to particle-hole symmetry, it is sufficient 
to consider < / < 1/2. The following distinctions are 
evident from Eqs. (3-5): 

(a) For low fillings (/ ^ 0,/i ^ -4t), a threshold 
coupling is required for d-wave pairing, while in the s* 
case, A^ ^ as F ^ due to a weak singularity at /i = 
—At. At 1/2-filling, however, due to a weak singularity 
at /i = in the d-wave case, A^ ^ as V ^ 0. For 5*, 
this singularity is not present, so, as A* 0, V/At 
7r^/8, i.e. a minimum coupling is needed for pairing. In 
contrast with A^ (V), A^{V) changes slope at /i = — 4t, 
and hence not smooth everywhere (though continuous). 

(b) For small momenta /c, the system exhibits the fol- 
lowing limiting behavior: (i) Ck < lJi{= — 4t); \uk\ 

1, |v/e| 0; this is the strong-coupling BEG limit. Here 
the ratio Vk/uk ~ A/c/2|/i| {k^ — /c^)/2|/i|, i.e. ana- 
lytic, (ii) e/e > /i(= — 4t); \uk\ 0, \vk\ 1; this is the 
weak-coupling BGS limit. Here Vk/uk ^/{kx — ky)^ i.e. 
non-analytic, (iii) Ck = fi{= — 4t); \uk\ 7^ 0, 7^ 0, 
when Ej^ 0. Then Vk/uk ^ {kx — ky)/{kx + ky)^ i.e. 
intermediate between (i) and (ii). For d-wave, the quasi- 
particle excitations in the BGS limit (ii) are "gapless" 
for some values of /c, while in the BEG limit (i), Ek 7^ 0, 
even for gaps with nodes [31]. 

Self-consistent numerical solutions of Eqs. (3-5) bear 
out the above features in detail, and reveal additional 
features. We scale /i, V, A by hopping parameter t. At 
a given filling /, both A^ and A^ increase with increas- 
ing V . While for d-wave it is easier to pair electrons 
at higher fillings, this is not necessarily the case for 5*- 
wave for the weaker couplings V/At < 1.5 and small gaps 

In Fig.l we show iiiV) for different fillings f. At a 
fixed /, in both the d- and 5*-wave cases, ji decreases 
with increasing coupling V, changing less rapidly for pro- 
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FIG. 1: Chemical potential /x versus coupling V at different 
fillings / for cZ-wave pairing. BEC pairs appear where IJ^{V) 
crosses the fi/2t = —2 line. The inset shows fi{V) for s- (dash- 
short dashed line), s*- (dashed line), and d-wave (solid line) 
pairing at / = 0.2. 

gressively larger /. However, as shown in the inset, for 
s*-wave, /i(V) exhibits a small "bump" for weaker cou- 
plings V/At < 1.5; for uniform 5- wave, the drop in fi 
with V is significantly more rapid. Crossover to the BEC 
regime is signaled by fJ^iV) going below the fi = —At line. 
As Fig.l shows, for d-wave, this develops at all fillings 
for some minimum coupling Vb/4:t . We note that as 
/ ^ 0, Vb/4:t 1.8; at 1/2-filling, this coupling tends 
to infinitely large values. For V > V^^ the system is 
conducive to BEC pairing, and for V < Vb^ the system 
exhibits BCS-like features. 

Fig. 2 shows the behavior of d-wave gaps as a func- 
tion of coupling V for different values of the chemical 
potential /i. The /i = —At curve represents the locus 
of Vb/At for different fillings (see Fig.l), and demarcates 
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FIG. 2: d-wave gap functions A/2t vs. nearest-neighbor cou- 
pling VjAt for different chemical potential /i. Inset: Results 
for s* case. \i = —At demarcates BEC and BCS regimes. 
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FIG. 3: (a), (b): 3D plots of d-wave fermion distribution 
functions vl. vs. kx — ky at filling / = 0.1, showing abrupt 
"jump" in vl. In BCS regime (a), /i = -3t, = 5.2t, 

V = 18. 7t, and in BEC regime (b), /i = -6t, = 0.6t, 

V = 5.2t. (c), (d): 3D plots of d-wave Vk/uk vs kx — ky 
for the same parameters as in ((a),(b)) respectively. In BCS 
regime (c) it can be seen to be no n- analytic; in BEC regime 
(d) it is analytic, (e), (f): The same as in (a), (b), but for 
s*-wave; the behavior is smooth. 

BEC and BCS -pair regimes. To the left is the /i > —At 
region wherein finite gaps of the BCS or intermediate 
BCS-BEC types exist. On a given constant-/i curve it 
may not be possible to have solutions for any arbitrary 
filling, but only those that satisfy Eqs. (3) and (4) self- 
consistently. The inset in Fig. 2 shows the corresponding 
A^ (y) curves for the s* case. There are interesting dif- 
ferences with the d-wave results in that the boundary 
(/i = —At) separating BEC/BCS regimes is not as clear- 
cut for weaker couplings V/At < 1.5 and smaller gaps 
A/2t < 0.5, however the /i < —At region lies to the right 
of the 11= —At curve, as in the d-wme case. 

Differences in the gap symmetry manifest in a strik- 
ing manner in the momentum distribution function, v^, 
and the ratio Vk/uk- For tZ-wave, for a given filling, in 
the weak-coupling BCS regime {V < Vb{f), ji > —At), 
vl exhibits a peak centered around the zone center (0,0), 
that becomes progressively narrow with decreasing fill- 
ing. Then at the crossover point at Vb{f) (/i = —At), 
vl abruptly goes to zero around (0,0), and shows a dras- 
tic redistribution along (0,zb7r), and (zLtt, 0) of BZ. The 
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abruptness is manifested in a "jump" in as the chemi- 
cal potential goes from just above the bottom of the band 
(/i > -At) to just below (/i < -4t), i.e. from BCS to BEC 
regime. A representative case is shown in Figs. 3a,3b. 
In marked contrast, for 5*-wave, (Figs. 3e, 3f), the zone 
center peak in v'j. decreases smoothly as one goes from the 
BCS regime to the BEC regime; only a slight redistribu- 
tion occurs at (zLtt, zLtt). This behavior is replicated at all 
fillings / < 1/2. As noted above in the small- /c limiting 
cases, our numerical calculations show (Fig 3c, 3d) that 
for d-wave, in the weak-coupling BCS regime, Vk/uk is 
non-analytic at ±:kx = ^ky; in the strong-coupling BEC 
regime, Vk/uk is analytic, vanishing along the zone diag- 
onals and peaking about (±7r,0), (0,±7r). In the s* case 
(not shown), v^/uk is analytic in both regimes. Thus 
we expect states with d- or 5* pairing gap symmetry to 
exhibit contrasting behavior at the BCS-BEC crossover, 
i.e. the unitarity limit. 

The Fourier transform of v'^ikx^ky), namely, py{x,y) 
reflect these differences. In the d-wave case, in marked 
contrast with its behavior in the BCS regime, py{x,y) 
is oscillatory in the BEC regime, and exhibits an inho- 
mogenous "checkerboard-type" pattern; see Fig 4(a,b). 
For the parameters of Fig. 3, the contrast ratio of the 
lowest density to the peak is roughly 50%, being most 
sensitive to the location of /i(F). The length scale is of 
the order of fractions of lattice spacing. py{x,y) is fairly 
uniform in the 5* case in both regimes. 

For correlated electrons, existent experimental tech- 
niques can reveal the proposed signatures to distinguish 
between states with d or extended-s (5*) symmetry; e.g. 
v'l, its Fourier transform, and quasiparticle energy could 
be deduced from ARPES [30], pair symmetries from 
quasiparticle tunneling or scanning tunneling microscopy 
(STM) In cold fermion systems, v^, and its Fourier 
transform, could be determined from time-of-flight [26] 
measurements, in which atoms are released from the lat- 
tice, and imaged at a later time. Combined with rf pair- 
ing gap spectroscopy [25l [33] (analog of tunneling), this 
could decipher the pairing symmetry. Momentum re- 
solved rf spectroscopy ^24, ^29^ , could provide information 

b 
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FIG. 4: (a) Fourier transform pv{x,y) of a typical d-wave 
fermion distribution function, v^. Here, filling / — 0.01, p — 
—4.2t (strong-coupling regime), gap A = .76t. (b) Projection 
of (a) to show contrast ratio of pv{x,y). 



on the quasiparticle energy Ej^ = (A^ + ^D^^^' thereby 
shedding light on {uk^Vk) and density of states. 

Here we have considered a homogeneous Fermi sys- 
tem, so as to appeal broadly to both condensed matter 
and cold atom physics. While we expect key aspects of 
our results to hold for cold fermion systems in the pres- 
ence of atomic traps, we comment on possible effects of 
the trap. But, flrst, we note that we consider a pair- 
ing Hamiltonian based on the extended Hubbard model, 
and not the simple Hubbard model, on which most dis- 
cussions in cold atom literature are based. We are also 
away from 1/2 fllling, and at relatively strong coupling, 
so that possible effects of spin density wave (SDW) and 
charge density wave (CDW) instabilities are expected to 
be suppressed. We have checked [14J that addition of 
next-near-neighbor hopping tends to stabilize the paired 
state, as well as lower the minimum near-neighbor inter- 
action necessary for a bound-state. 

Away from 1/2 fllling (lower flllings) for a 2D lattice 
in presence of trap (within local density approximation 
(LDA)), the multitude of quantum phases obtained [34] - 
|36] near 1/2 fllling may not exist (34[ [35] , leaving a metal- 
lic state from the trap center to the edge. Thus, for this 
range of flllings, our results are not expected to be sub- 
ject to competing signatures from other possible phases 
in rf spectra. Closer to 1/2 fllling, within LDA, there 
may be other phases, e.g. pairing in different shells be- 
having as inflnite system in each shell [3^. Then, tomo- 
graphic (spatially resolved) rf spectroscopy ^ , using in 
situ phase contrast imaging technique, would be able to 
probe each shell region, within each of which our results 
should hold. 

Recent work [37 have pointed out possible effects of 
trap inhomogeneity on the Hartree term, and conse- 
quently on rf spectra. In this work, Hartree effect on 
self-energy is obtained at flrst-order, so the role of higher- 
order terms towards a convergent result is not clear at 
this point. It may also be interesting to examine Hartree 
effect for the extended Hubbard model on which our cal- 
culation is based; this is outside the scope of current 
work. 

Our calculation is at T= 0, and any Berezinski- 
Kosterlitz-Thouless transition [38 would only be re- 
vealed in a flnite-T calculation. Also, like others, we 
take our the system to be not strictly 2D, but quasi-2D, 
with the assumption that a weak link along the 3rd di- 
rection stabilizes phase transitions like superfluidity. Our 
T=0 consideration does not lend itself to calculations 
of critical temperature, Tg. However, based on other 
work [lot I39j , for d-wave pairing we estimate that in elec- 
tron systems, such as the cuprates, Tc/Tp ~ 0.015 — 0.03, 
giving a Tc - IbK - 30i^ for Tp ^ lO^iC; in cold 
fermions, Tc/Tp ~ 0.01, giving a Tc ~ 30ni^ for a 
Tp ^ 3/ii^. Thus for d-wave Tc measurements in cold 
fermions, realizing temperatures below the currently at- 
tainable T/Tp 2^ 0.05 is needed. However, these are 
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lower bound estimates, and we expect our proposed sig- 
natures to persist to higher temperatures. Recent sug- 
gestions [40 of novel cooling methods are encouraging. 

Our calculations, in the spirit of BCS and BEC-BCS 
crossover theories, consider d- and extended-s (s*) wave 
pairing symmetries, independent of pairing mechanisms. 
Though mean-field in nature, we expect the calculated 
signatures of unconventional pairing symmetries to hold 
in calculations beyond mean fields. A recent Monte Carlo 
work [2T| on fermions on 2D optical lattice, though at 1/2 
filling, goes beyond mean field and contain substantial 
discussions regarding possible phases. 
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